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ABSTRACT 

We have identified CFBDSIR J1458+1013 as a 0.11" (2.6 AU) physical binary us- 
ing Keck laser guide star adaptive optics imaging and have measured a distance of 
23.1 ±2.4 pc to the system based on near-IR parallax data from CFHT. The integrated- 
light near-IR spectrum indicates a spectral type of T9.5, and model atmospheres suggest 
a slightly higher temperature and surface gravity than the T10 dwarf UGPS J0722— 05. 
Thus, CFBDSIR J1458+1013AB is the coolest brown dwarf binary found to date. Its 
secondary component has an absolute iif-band magnitude that is 1.9±0.3 mag fainter 
than UGPS J0722— 05, giving an inferred spectral type of >T10. The secondary's bolo- 
metric luminosity of ~ 2 x 10~ 7 L Q makes it the least luminous known brown dwarf by 
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a factor of 4-5. By comparing to evolutionary models and T9-T10 objects, we estimate 
a temperature of 370±40 K and a mass of 6-15 M Jup for CFBDSIR J1458+1013B. At 
such extremes, atmospheric models predict the onset of novel photospheric processes, 
namely the appearance of water clouds and the removal of strong alkali lines, but their 
impact on the emergent spectrum is highly uncertain. Our photometry shows that 
strong CH4 absorption persists at -ff-band; the J — K color is bluer than the latest 
known T dwarfs but not as blue as predicted by current models; and the J — H color 
delineates a possible inflection in the blueward trend for the latest T dwarfs. Given 
its low luminosity, atypical colors and cold temperature, CFBDSIR J1458+1013B is a 
promising candidate for the hypothesized Y spectral class. However, regardless of its 
ultimate classification, CFBDSIR J1458+1013AB provides a new benchmark for mea- 
suring the properties of brown dwarfs and gas-giant planets, testing substellar models, 
and constraining the low-mass limit for star formation. 

Subject headings: binaries: general, close — stars: brown dwarfs — infrared: stars - 
techniques: high angular resolution 



Introduction 



Over the last 15 years, stellar and exoplanetary studies have been propelled through the study 
of brown dwarfs. Residing at the extremes of low mass, luminosity and temperature, brown dwarfs 
serve as laboratories for understanding gas-giant extrasolar planets as well as the faint end of the 
star formation process. The coolest known brown dwarfs, the T dwarfs, have temperatures (~600- 
1400 K), bolometric luminosities (~ 10~ 4 to 10~ 6 Lq), molecule-dominated (H2O and CH4) IR 
spectra, and atmospheric processes (e.g., non-equilibrium chemistry) that are more akin to Jupiter 
than any star. This physical similarity drives the ongoing search for even more extreme objects, to 
close the gap of «400 K in temperature and factor of «1000 in L\ m \ between the coolest T dwarfs 
and the planet Jupiter. Along the way, one question to be answered is whether su ch new discoveries 



will h ave sufficiently distinct spectra to invoke the proposed "Y" spectral type (jKirkpatrick et al 
19991 ). 



A well-established strategy for finding the coolest objects is to search for them as compan- 
ions to other objects. This approach has been quite s uccessful: the first L and T dwarfs were 
found as companions around hotter objects, GD 165B (jBecklin fc Zuckermanlll988l ) and Gl 229B 
(jNakaiima et al.lll995l ). respecti vely. To date, some of the coolest T dwarfs are found as compan- 
ions to stars: Gl 57 D (T7.5; bureasser et aP l200oh. HP 3 651 B fT7.5: llVlugrauer et alJ bood: 
Luhman et alJl2007h. UL AS J141 6+13 B (T7.5: ISchold boiO^l: IBurningham et alJboioh. GJ 758B 



(«T8-T9; iThalmann et al 



and Ross 458C (T8.5; 



2009 



Janson et al 



Goldman et al 



2010: 



Scholz 



201 ll) Wo lf 940B (T8.5; IBurningham et al.l 120091 ) 



2010b! ) . 



Field L and T dwarfs themselves have long been recognized as promising search targets for even 
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cooler objects. To date, six second ary components later than T5 have been identified in binaries 
(see compilation in iLiu et al.l [201C]) , one of which has a short enough orbi tal period to y ield a 



dynamical mass, the T5.0+T5.5 substellar binary 2MASS J1534-2952AB (|Liu et al.ll200a )FI In 



addition, discovery of ultracool objects in binaries opens the door to direct mass measurements 
through astrometric and spectroscopic monitoring of their orbits. 

As part of an ongoing program to study the multiplicity and physical properties of brown dwarfs 
using laser guide star adaptive optics (LGS AO), we present here the discovery of the binarity of 
CFBDSIR J1458+1013. This object was identified in deep z and J-band data by iDelorme et al 



(|2010l ) and assigned a preliminary spectral type of T8.5 based on its H-band spectrum. We have 
scrutinized this object using a combination of high resolution imaging, near-IR astrometry, and 
near-IR spectroscopy. These data indicate an integrated-light spectral type of T9.5, wh ich is 
later than all known field objects except the T10 dwarf UGPS J0722— 05 (ILucas et al.l 12010 ). The 
secondary component of CFBDSIR J1458+1013AB is ~1.5-2 mag fainter in absolute magnitude 
and «150 K cooler in temperature than any known brown dwarf, making it a promising candidate 
for the first member of the Y spectral class. 



2. Observations 
2.1. Keck LGS AO Imaging and Photometry 

We imaged CFBDSIR J1458+1013 on UT 22 May 2010 and 08 July 2010 using the sodium 



LGS AO system of th e 10-meter Keck II Telescope on Mauna Kea, Hawaii (|Wizinowich et ajj|2006| ; 



van Dam et al.ll2006l ). We used the facility IR camera NIRC2 with its narrow field-of-view camera, 
which produces a 10.2" x 10.2" field of view. Conditions were photometric for both runs with 
excellent seeing: 0.5-0.6" and 0.3-0.5" FWHM in May and July, respectively, as measured by a 
differential image motion monitor located at the nearby Canada-France-Hawaii Telescope. The 
LGS provided the wavefront reference source for AO correction, with the tip-tilt moti on measured 



conte mporaneously from the R = 15.5 mag field star USNO-B1.0 1002-0236935 (jMonet et al 



20031 ) located 55" away from CFBDSIR J1458+1013. The LGS brightness, as measured by the flux 
incident on the AO wavefront sensor, was equivalent to a V ~ 9.6 and 10.3 mag star in May and 
July, respectively. 

At each epoch, we obtained a series of dithered images by offsetting the telescope by a few 
arcseconds between each 1-2 images. The sodium laser beam was pointed at the center of the NIRC2 
field-of-view for all observations. In May 2010, we successfully resolved CFBDSIR J1458+1013 into 
a tight binary using the moderate-band CH^s filter, which has a central wavelength of 1.592 /jm 
and a width of 0.126 /jm. In July 2010, we obtained images with the J (1.25 /xm), H (1.64 /um), 



The T8 dwarf 2MASS J0939— 24 ha s also been suggeste d to be an unresolved binary based on its absolute 
spectrophotometry ( Burgasser et al. 20081 ; Leggett et al. 20091 ) . 
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and K (2.20 tim) filters from the Mauna Kea Observatories filter consortium (pinions Tokunaga 
20021 : iTokunaea et alJliooj ). 



The images were analyzed in a similar fashion as our previous LGS papers (e.g. lLiu et al.l l2006. 



2008). The raw images were reduced using standard methods of flat-fielding and sky-subtraction. 
After rejecting individual images of poor quality, we registered and stacked the data into a final 
mosaic. The binary's flux ratios and relative astrometry were derived by fitting an analytic model 
of the point spread function (PSF) as the sum of two elliptical gaussians to the images of the binary 
in the mosaics. We did not correct the relative astrometry for instrumental optical distortion, as the 
size of the effect as predicted from a distortion solution by B. Cameron (priv. comm.) is insignificant 
compared to our measurements uncertainties. We adopted a pixel scale of 9.963 ± 0.005 ma s/pixel 
and an orientation for the detector's +y axis of +0.13 ± 0.07° for NIRC2 (jdiez et alJbood ). 



To validate our measurements, we created myriad artificial binary stars from the science 
data themselves; given the large flux ratio of the binary, a sc aled and shifted ve rsion of the 
primary star provides an excellent simulation of the data (e.g. iDupuv et al.l l2009bl ). For each 
filter, our fitting code was applied to artificial binaries with similar separations and flux ratios as 
CFBDSIR J1458+1013AB over a range of position angles (PAs). The RMS scatter of the sim- 
ulations were adopted as our measurement errors; the results also showed that any systematic 
offsets are small. The exceptions to this were the J and i?-band separation measurements, where 
the average systematic offset from the simulations was slightly larger than the RMS scatter, and 
for these datasets we adopted the offset as the uncertainty. Table [Q presents our final Keck LGS 
measurements, and Figure Q] shows our Keck LGS images. 



We use our measured flux ratios and the published JHK$ photometry from iDelorme et al 



(|20ld ) to derive resolved IR colors and magnitudes for CFBDSIR J1458+1013AB on the MKO 
system. Delorme et al. report data using the -f^-band filter while our Keck LGS imaging used 
the i^-band filter; to compute resolved i^-band photometry, we synthesize an integrated-light 
color of (K - Ks) = -0.17 ± 0.02 mag from the published spectra of ULAS J1238+09 (T8.5), 
Wolf 940B (T8.5), ULAS J0034-00 (T9), CFBDS J0059-01 (T9), an d ULAS J1335+11 (T9) 



(jBurningham et al .1 120081 . 



2009 



Warren et al. 



2007 



Delorme et al.l l2008al ) . This is consistent with 



the color synthesized from the CFBDSIR J1458+1013AB spectrum itself (§ I2~3]) . but the signal- 
to-noise (S/N) in the iT-band spectrum is quite low, so we use the average from the known T8.5 
and T9 objects. To obtain the resolved CH^s photometry from the published ff-band data, we 
synthesizes a color [CH^s — H) = —0.47 ± 0.03 mag directly from the integrated-light spectrum of 
CFBDSIR J1458+1013AB. Table H presents the resolved photometry. 



2.2. CFHT Infrared Astrometry 



To determine the distance to CFBDSIR J1458+1013AB, we acquired seven epochs of J-band 
imaging in 2009 and 2010 spanning 1.07 years using the Canada- France-Hawaii 3.6-meter Telescope 
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on M auna Kea, Hawai'i. We used the facility wide-field near-IR camera WIRCam (jPuget et al 
2004 ). which has four 2048x2048 Hawaii-2RG detectors, each with an approximate field-of-view of 
10' x 10'. At each epoch, we placed the target near the center of the northeast detector. The data 
were taken with slightly different observing procedures as they originated from two independent 
programs: at 3 of the epochs, images of 60-s x 6 coadds were obtained at 10 dither positions, 
and at the other 4 epochs, single 60-s images were obtained at 11 dither positions. The median 
FWHM of the data was / .'73 (min/max of / /55/0 / . / 80). The target was always observed near transit, 
with a maximum deviation in airmass of 0.03. This reduces any systematic offsets in the target's 
astrometry from the differential chromatic refraction between the late-T dwarf science target and 
the background stars to < 1 mas (Dupuy & Liu, in preparation), well below the measurement noise. 



Our WIRCam astrometry pipeline is described in detail by lDupuyl (120101 ) and Dupuy & Liu (in 
prep) and successfully generates par allaxes that are in accord with previously published results for 
L and T dwarfs. We use SExtractor (jBertin &; Arnoutslll996l ) to generate an astrometric catalog of 
all the sources in every individual image. We then cross-match detections across dithered images 
and use the mean and rms of each source to determine its final position and error at that epoch; 
we exclude any source that has fewer than 7 detections or a S/N less than 7. We then cross-match 
these positions across different epochs, allowing for a full linear transformation in the astrometry 
to account for scale and orientation changes in the instrument. (WIRCam is installed and removed 
from the telescope for each observing run.) After an initial iteration, we remove any stars with 
large proper motions (> 50 mas/yr) or parallaxes (> 3a) from the epoch-matching fit. Finally, 
we absolutely calibrate our astrometry by matching to 2MASS to determine the overall scale, 
orientation, and shear. For this target, 8 of our 72 reference stars were well-detected in 2MASS 
and had sufficiently low proper motion (< 30 mas/yr) to be used for absolute calibration. 

We then fit for the proper motion and parallax of every source in the field. We measured 
a proper motion of 432 ± 6 mas/yr at a PA of 154.2 ± 0.7 deg for CFBDSIR J1458+1013, in 
good agreement with the lDelorme et al .1 (|2010l ) measurement of 444 ± 16 mas/yr and 157.5 ± 2.1° 
from seeing-limited images. CFBDSIR J1458+1013 was the only object we measured in the field 
with a significant parallax, 42.4 ± 4.5 mas relative to the reference stars. (To our knowledge, 
this is the faintest objec t outside the solar system with an optical or IR parallax measurement, 
cf., iMarocco et al.ll2010l .) To a ccount for the me an parallax of our reference stars, we used the 
Besangon model of the Galaxy (jRobin et ajj|20fl3j ) to simulate the distance distribution of sources 
in the direction of CFBDSIR J1458+1013 within the magnitude range of our images (15.1 < J < 
20.1 mag). We find a mean parallax of 0.9 ± 0.1 mas from the Galaxy model, with the error 
accounting for the sampling variance due to the finite number of reference stars. This results in an 
absolute parallax of 43.3 ± 4.5 mas for CFBDSIR J1458+1013 ( 2 3.1 ± 2.4 pc), in good agreement 
with the original 23 pc photometric estimate by iDelorme et al.l (|2010l ). The resulting tangential 
velocity (47±5 km s _1 ) is higher than the median for field T dwarfs, but not extraordinary give n 
the observed velocity dispersion of the field objects (e.g. IVrba et al.1 120041 : iFaherty et al.ll2009i ). 
Table [2] contains our astrometry results. 
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2.3. VLT/X-Shooter Spectroscopy 

We obtained integrated-light near-IR spectra of CFBDSIR J1458+1013AB with the X-Shooter 
spectrograph on VLT-UT2 (Program 085.C-0805). The observations were carried out in four 1-hour 
observing blocks (hereafter OB) from May 5 to July 9, 2010, achieving a total exposure time on 
target of 11700 sec, split into 8 A-B nods on the slit of 2 x 732 sec each. The slit width was 1.2", 
and the seeing varied between 0.8-1.0". We chose a large slit size, because direct acquisition of the 
target was impossible and we had to acquire using a blind offset from a reference source. 



The spectra for each OB were reduced using the standard ESO X-Shooter pipeline (jModigliani et al 



2010h . which produced a 2-dimensional, curvature-corrected spectrum of the NIR arm of X-Shooter 
(1.0-2.5 //m). No signal was detected in the UVB and visible arm of the instrument. The trace 
was extracted using our own IDL procedures, using gaussian boxes in the spatial dimension at each 
wavelength. A similar gaussian extraction box was used at about 5 times the FWHM of the trace 
to obtain the noise spectrum of the sky The resulting 1-dimensional spectrum from each science 
target OB was then divided by the spectra of telluric stars observed just before or after, reduced 
and extracted using the same pipelines. 

The resulting 4 spectra from all the OBs were then median combined to produce the final 
science spectrum. Since the resulting spectrum (with R = A/ A A ~ 4300) has low S/N and most 
analyses of late-T dwarfs have used much lower resolution spectra, we smoothed the spectra using a 
weighted average over 50 pixels in the wavelength dimension (corresponding to 3 and 5 nm at J-band 
and JC-band, respectively), weighting by the inverse variance of the extracted sky spectrum. This 
approach uses the full spectral resolution of X-Shooter to give a much lower weight to wavelengths 
affected by telluric emission lines, thus improving the S/N with respect to simple unweighted 
binning using the average or median. 

Given the large wavelength range covered by X-Shooter, seeing variation in the spectral direc- 
tion can led to wavelength-dependent slit losses. Th erefore, we flux - calibr ated the spectrum using 



JHK S photometry of CFBDSIR J1458+1013 from iDelorme et al.l (|2010h and F-b and data from 



the UK IDSS Data Release 8 (Y = 20.58 ±0.21 mag, on the Vega system). Following IDelorme et al 



(|2008bl ). we applied small scaling factors to each of the broadband wavelength ranges so that the 
magnitudes derived from the CFBDSIR J1458+1013 spectrum matched the photometry. The final 
spectrum is shown in Figure [H 



3. Results 



3.1. Proof of Companionship 



Examination of the Digitized Sky Survey images shows no plausible optical counterpart if 
co mponent B wer e a ba ckground object. Similarly, the z-band imaging from July 2004 presented 
Delorme et al.l (|2010l ) shows no star to a depth of z ~24 mag, 2-3 mag deeper than expected 



m 
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for a background M star if it had the near-IR magnitude of component B. In fact, the high proper 
motion of the source means that our two Keck epochs separated by 1.5 months are sufficient to 
confirm the two components are physically associated. (For the July 2010 epoch, we choose the 
Keck ii-band astrometry, which has the smallest errors.) Figure [2] shows the large expected motion 
for component B if it were an object at infinite distance, given our measured proper motion and 
parallax for CFBDSIR J1458+1013. (The proper motion dominates over the parallax.) In contrast, 
the relative positions of the two components do not change between our observations. The length 
of the vector between the observed and predicted change is non-zero at 9.1a. Thus, we conclude 
the two components are physically associated. 



3.2. Integrated-Light Spectrum 

3.2.1. Comparison with known late-T dwarfs 

The integrated-light spectrum of CFBDSIR J1458+1013AB is very similar to other recent very 
late-T dwarf discoveries (Figured]). As seen by its narrower continuum peak, its J-band spectrum 
shows slightly stronger absorption than the T9 dwarfs CFBDS J0059-01, ULAS J0034-00, and 
ULAS J1335+11 and is very close to the T10 dwarf UGPS J0722-05@ (We caution that spectral 
types later than T8 are an active topic of investigation, without a classification system established 
yet. For this discussion, we simply adopt the published types for these objects, which may ultimately 
change once a system is defined.) In the H band however, CFBDSIR J1458+1013 is more similar 
to the warmer objects, especially on the blue side of the peak. The absorption on the red side of 
the i7-band is the strongest observed in any known brown dwarf (see the CH4-H index in Table H]) . 
This might be due to enhanced collision induced absorption (CIA) of H2 starting to suppress the H- 
band peak; this notion fits well with the strongly depressed X-band flux of CFBDSIR J1458+1013 
compared to the other late-T dwarfs, where the strongest continuum absorber is CIA H2. 

The suppression could be due to e ither strong gra vity or low metallicity, and d istinguishing be- 



tween t he t wo can be chal l enging (e.g. iLiu et al.l 120071 ). Using model atmospheres, iBurgasser et al 



(|2006al ) and lLeggett et al.l (|2007l ) show that the shape of the Y-band continuum is more metallicity- 
dependent than gravity-dependent (though their two sets of models show different amplitudes for 
these effects). The main absorbers on the blue side of the Y-band peak are the pressure broadened 
wings of the K I and Na I fundamental doublets, which naturally decrease with lower metallicity; 
more metal-rich models also have the Y-band peak at redder wavelengths and the overall Y-band 



Lucas et al.l l |20ich report the presence of an unidentified, narrow absorption feature at 1.275 /im in the spec- 



trum of UGPS J0722— 05. They suggest it arises from H2O or HF. Figure U suggests it may also be present in 
CFBDSIR J1458+1013AB. However, while the detection is significant in terms of S/N, it appears much broader 
than in the UGPS J0722— 05 spectrum. Also, Lucas et al. reports a second absorption feature at 1.282 /im, which 
coincides with a local peak in the CFBDSIR J1458+1013AB spectrum. Thus our spectra do not confirm either of 
these proposed features. 
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flux being fainter. The observed Y-band spectra of CFBDSIR J1458+1013AB does not show any 
strong differences with UGPS J0722-05 and CFBDS J0059-01 (Figure HJ), though the peak might 
be slightly redder than the other two objects, perhaps hinting at a greater metallicity (though 
this would be discrepant with the highly suppressed -ff-band). Overall, we favor the higher gravity 
explanation, though a better spectrum at shorter wavelengths would help to investigate this further. 

Taken altogether, the spectrum of CFBDSIR J1458+1013AB shows that this object is slightly 
earlier than UGPS J0722— 05 and possibly higher gravity, but later than the known T9 dwarfs. 
Thus, we provisionally assign a spectral type of T9.5, though of course it may turn out to be earlier 
or later type, depending on the future development of classification for such cool objects. 



3.2.2. Spectral index analysis with solar-metallicity models 



To quantify the absorption features of CH 4 , H oO, CIA H?, and poss i bly NH3, we calcu- 



lated spectral indices previously used for T dwarfs (jBurgasser et al.l l2006bl : IWarren et al. 



2007 



Delorme et all l2008al ). The measurements (Table 0$ indicate that CFBDSIR J1458+1013A has 
some of the strongest absorption features seen in the near-IR spectra of late T dwarfs, exceeded 
by only UGPS J0722— 05 for the Wj and NH3-H indices. Given the «2 magnitudes flux difference 
between component A and B, it is reasonable to consider the integrated spectrum is dominated by 
the primary (as discussed below). 



We use the solar-abundance BT-Settl-2009 models (jAllard fc Freytag! l2010l ) to estimate T e g 
and log(p) for com ponent A from the index measureme nts. We compensate for the systematic errors 
of the models (e.g. iLeggett et al.ll2007l : iLiu et al.ll2008l ) by scaling their index pr edictions to match 
those measured for the T8.5 benchmark Wolf 940B (jBurningham et al.l 120091 ). Its high quality 
spectrophotometry, parallax, and age estimate from its primary star enable a pr ecise estimate of its 
prop erties usi ng evolutionary model s (T e g = 575 K; log(g) = 4.75, [M/H]=0 from lBurningham et al. 
20091 : see also ILeggett et al.ll2010al ). With this approach, the temperatures and gravities for the 



>T8 dwarfs derived from an index-in dex plot broadly agree with the values derived from the 



analysis of their full near-IR spectra (jBurgasser et al.l l2006al : IWarren et al.1 120071 : iDelorme et al 
2008al : iBurningham et al.l boosl : Licas et alJboid ) . 



Since the systematic errors in the models in this low-temperature regime are likely much larger 
than the measurement errors, we use the model grid only for relative estimates of temperature 
and gravity differences among the late-T dwarfs, not absolute measurements. Figure [5] indicates 
that CFBDSIR J1458+1013A is as cool as UGPS J0722-05 and >50 K or more cooler than the 
other >T8 dwarfs. (The caveat is that the model analysis assumes a single metallicity for all 
the objects, consistent with the similarity in Y-band shapes discussed above.) The plot also puts 
CFBDSIR J1458+1013A at a higher (>0.5 dex) surface gravity than the other objects, as expected 
from the highly suppressed J^-band flux. The higher gravity suggests that CFBDSIR J1458+1013A 
is older and higher mass than UGPS J0722— 05, consistent with their relative tangential velocities 
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(47±5 km s 1 and 19±3 km s , respectively). Of course, tangential velocities are only a statistical 
proxy for age so this comparison should be treated accordingly. 

To assess the impact of component B on the integrated-light spectrum, we used the BT-Scttl- 
2009 model atmospheres to create composite spectra. We chose models with 600-700 K for the 
primary and 400-500 K for the secondary, matching pairs of models with 200-250 K temperature 
differences and the same surface gravity. The pairs of spectra were scaled to agree with the observed 
J-band or U-band flux ratio and summed to form artificial composites. We then added noise 
based on the measured S/N of the VLT spectrum and measured the same indices. We found that 
the blended composites had very similar index values as the input primary spectrum alone, with 
systematic effects at the 5-20% level, typically 10%. Such amplitudes are not sufficient to affect 
the interpretation of Figure [5j 



3.2.3. Model atmosphere fitting 



To further place CFBDSIR J1458+1013AB in context, we perform a homogeneous comparison 
of the near-IR spectra of very late-T dwarfs to low-temperature atmospheric models. We restrict 
our analysis to >T8 dwarfs with parallaxes, excluding the aforementioned possible binary T8 dwarf 
2MASS J0939-24, but also include two benchmark T7.5 dwarfs (Table We use the published 
J-band photometry to establish the absolute flux calibration of the spectra, again assuming the 
primary dominates the integrated-light spectrum. For CFBDSIR J1458+1013AB, we again as- 
sumed that the integrated-light spectrum is dominated by component A and normalized it using 
the resolved J-band photometry of component A. 



2001 



We fit three grids of atmospheric models: (1) the Ames-Cond models (jAllard et al 
spanning T eS = 400 - 1200 K (Ar pff = 100 K) and log(s) = 4.0 - 6.0 (Alog(s) = 0.5); (2) the 



Tucson models (jBurrows et al.ll2003l . hereafter BSL03), which are irregularly gridded between T c g « 
150 — 850 K and log(g) ~ 3.3 — 4.8 (cgs); and (3) the aforementioned BT-Settl-2009 models spanning 
400-1200 K and log(g) =4.0-6.0 dex (AT eff = 50 K above 600 K and 100 K below 600 K). All 
grids are for solar metallicity compositions and chemical equilibrium. The models were gaussian 
smoothed to the resolution of the data and interpolated onto the wavelength grid of the observed 
spectrum. 

The atmospheric models have known systematic errors, particularly in the l"-band and the 
red side of the /J-band as a result of imperfect modeling of the pressure-broadened K I doublet 
in the far-red and the incomplete list of CH4 lines, respectively. We therefore restrict our fits to 
the 1.15-1.35 /im, 1.45-1.6 /jm, and 1.95-2.25 /im spectral regions. The model emergent spectra 
ar e scaled to the flux calibrated spectra by minimizing the x 2 statistic (e.g., Equation 5 and 6 
of iBowler et al.ll2009l ). For observed spectra without measurement errors (C FBDS J0059— 01 and 



UGPS J0722— 05) we perform a least-squares fit to the data (Equation 1 of IStephens et al 
with Wi=l, and the scaling factor is solved for by equating the derivative of G\ with zero). 



2009 
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The best-fitting models are summarized in Table [5] and plotted in Figure [6j The Ames-Cond 
models yield effective temperatures between 700-800 K and a gravity of 4.5 dex for almost all 
objects; the temperatures for the latest-type objects are ~100-200 K hotter than analyses from 
their discovery papers (which use variants of the BT-Settl models). The BSL03 models typically 
produce much cooler temperatures (420-800 K) and somewhat higher gravities (4.6-4.8 dex), with 
the results for the earlier-type objects often at the edge of the model grid, indicating a broader set 
of models is probably needed. The BT-Settl-2009 model results are usually intermediate between 
those from the other two models. The fits show that CFBDS J0059-01, CFBDSIR J1458+1013AB, 
and UGPS J0722— 05 are the three coldest objects in the sample, with the Ames-Cond models 
indicating all have comparable temperatures while the BSL03 and BT-Settl-2009 models indicate 
that CFBDSIR J1458+1013AB is the warmest. (Note that in terms of the absolute values derived 
from the fitting, the analysis in the Appendix on the resulting spectroscopic distances suggests that 
the temperatures may be somewhat overestimated.) 

Again, to assess the effect of component B on the analysis of the integrated-light spectrum, 
we took the artificial noisy composite spectra described in § 13.2.21 and applied the same model 
atmosphere fitting procedure. We found that the T e g's determined for the composite spectra were 
at worst 50 K colder than those of the input primary spectrum alone, and often the same. The 
fitted surface gravities were the same for the composites and the primary spectrum. Thus, the 
blended light from the much fainter secondary appears to have little consequence for our analysis, 
as expected. 

We also examine another method to derive physical properties from atmospheric models. Al- 
though the distance estimates using the scaling factor (R/d) 2 may differ substantially from the true 
distances, the known distance can be used in combination with the scaling factor, a grid of radii, 
and a likely age range to infer a temperature and mass range for the object. This is effectively using 
the scaling factor as a proxy for temperature, and it has the benefit of relying on the absolute level 
of the model emergent flux density rather than the detailed shape of the synthetic spectrum, which 
may have systematic errors. Using the results from Tabled we compute the range of effective tem- 
peratures and masses for the three sets of models. The results are shown in Table El assuming an age 



range of 1-10 Gyrs, except for Ross 458C where we assume an age of 0.1-1.0 Gyr (j Goldman et al 



2010). The BSL03 models are limited at high gravities by an upper mass cutoff of 25 Mj up so they 



only provide lower limits on T c g- and mass. The resulting r ange of temperatures and masses is sim 



ilar to those estimated in previous work (jLiu et al.l 120071 ; iBurningham et al.l 120081 ; iDelorme et al 



2008al ) This analysis indicates that CFBDSIR J1458+1013A has a similar effective temperature 
(540-660 K) and mass (12-42 M Jup ) to the T9 dwarfs ULAS J0034-00 and ULAS J1335+11, with 
CFBDS J0059— 01 and UGPS J0722— 05 being perhaps slightly cooler, in accord with the relative 
ordering of their absolute magnitudes. 
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3.3. Resolved Colors and Magnitudes 



In accord with the integrated-light spectrum, the near-IR colors of component A are similar 
to those of the known T8-T10 dwarfs, with its (H — K) color of —0.45 ± 0.26 mag being somewhat 
bluer than most objects, though the uncertainty is large. As discussed above, the bluer color is 
likely due to a higher surface gravity (i.e., an older age) compared to most other late-T dwarfs. 

Our near-IR parallax for the system provides the absolute magnitudes for component B and re- 
veals its extreme nature. Compared to the T10 dwarf UGPS J0722-05, CFBDSIR J1458+1013B is 
fainter by 1.4±0.4, 1.9±0.3, and 2.0±0.4 mag in its JHK absolute magnitudes, respectively. Being 
significantly fainter than previously known brown dwarfs (Figure[7|), the properties of CFBDSIR J1458+1013B 
provide a first look at the empirical trends beyond spectral type T10, especi ally enlightening give n 
the uncertainties in model atmosphere predictions for such cold objects (e.g. iBurrows et al.ll2003l ). 
The very blue (CH^s — H) color indicates strong CH4 absorption persists at f/-band, comparable 
to the known T8- T9 dwarfs which have a (CH^s — H) ~ —0.5 to —0.6 mag (e.g., Figure 2 of 
Tinnev et aDl2005l and Section 3.2 of lLiu et aliboosh . 



The near-IR colors of CFBDSIR J1458+1013B are very blue (Figure ED, with its H - K color 
is similar and the J — K color slightly bluer than previously known T dwarfs (except for th e T6p 
dwarf 2MASS 0937+29, which is suspected to have low metallicity; iBurgasser et al.ll2002l ). For 
i ts absolute magnit ude, its J — K color is much redder than those predicted by the Lyon/Cond 
(jBaraffe et al.ll2003l ) and ISaumon &, Marleyl (|2008l . hereafter SM08) models, and to a lesser degree 
the BSL03 models as well. The most intriguing behavior is seen in J — H, where the blueward 
trend with spectral type among the T dwarfs appears to reach an inflection point with the T8-T10 
dwarfs, and then CFBDSIR J1458+1013B continues to bluer colors. 



3.4. Physical Properties from Evolutionary Models 



How much cooler is CFBDSIR J1458+1013B than previously known late-T dwarfs? Using its 
absolute magnitudes as a simple proxy for L^y (i.e., assuming the same bolometric corrections as 
other T9-T10 objects), its L^oi is a factor of ~4-6x smaller. Old brown dwarfs have essentially 
constant radii, so this would correspond to an object that is a factor of ~1.5 cooler. The T9-T10 
objects ha ye estimated T P ff of 500-600 K from model atmosphere fitting an d their bolometric lu- 



minosties (jLeggett et al.ll2009l . 



2010a 



Marocco et al. 



2010 



Lucas et al.ll2010l ): simple scaling would 



then make CFBDSIR J1458+1013B £3350-400 K. 

With Lboi estimates for the individual components and an assumed age, evolut ionary models 



can f ully determine the remaining physical properties of the two components (e.g. ISaumon et al 



2000). This approach is desireable as predictions of L^ \ are less dependent on the assumed model 
atmo spheres than predictions of the absolute magnitudes (jChabrier et al.ll2000bl : ISaumon &: Marley 
20081 ). This approach has also been widely used for temperature determinations for field L and 
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T dwarfs (e.g. iGolimowski et al.l 12004 ). Of course, such estimates should be treated with caution 
since CFBDSIR J1458+1013B resides in a low-temperature, low-luminosity reg ime where the mod- 



els have not been tested against any observations. (In the case of L dwarfs, iDupuv et al.l 12009a 



have found that for a given luminosity and assumed age, the evolutionary models may overestimate 
the masses.) 

To estimate Lbob we adopt J-band bolometric corrections from the coolest T dwarfs with 
parallaxes: ULAS J0034-00 (T9; BCj = 1.96± 011 mag), CFBDS J0 059-01 (T9; 1.87±0.21 mag), 
and ULAS J1335+11 ( T9; 1.82±0.0 9 m&g) from lMarocco et al.l (|20ld ). and UGPS J0722-05 (T10; 
1.36 ± 0.21 mag) from iLucas et all (|2010l ). The unweighted average is 1.75 ± 0.27 mag, and we 
apply this to both components of CFBDSIR J1458+1013AB. The resulting values for L^ \/Lq are 
(1.1 ± 0.4) x 10 -6 and (2.0 ± 0.9) x 10~ 7 for components A and B, respectively. In comparison, 
CFBDS J0059-0114 and UGPS J0722-05 have L ho i/L Q = (7.2 ±1.2) x 10" 7 and (9.2 ±3.1) x 10~ 7 , 
respectively, making CFBDSIR J1458+1013B a factor of ~4-5 less luminous than these two objects. 

Table [3] presents the derived physical quantities for ages of 1.0 and 5.0 Gyr. The Lboi un- 
certainty is propagated into the calculations in a Monte Carlo fashion. The inferred temperature 
of compo nent B is exce p tiona lly low, 360-380 K from the Lyon/COND models and 350-380 K 
from the Burrows et al. (jl997l ) ones. Including the random uncertainties from the Monte Carlo 
results, we adopt a final T e g estimate of 370±40 K. The model-derived mass is 6-14 (7-17) for 
the Lyon/COND (Burrows) models, possibly placing the object below the notional 13 Mj up di- 
viding line between planets and brown dwarfs based on the deuterium-burning limit. Overall, the 
physical properties of CFBDSIR J1458±1013B overlap those of the gas-giant planets found from 
radial velocity surveys. In fact, both components of CFBDSIR J1458±1013AB would reside in the 
planetary-mass regime and have retained their initial deuterium abundance for ages of <1 Gyr. 

To gauge the systematic effects, we also compute physical properties using the BSL03 models, 
which provide predictions for absolute magnitudes. (The Lyon/COND models also provide absolute 
magnitudes, but these disagree more strongly with the color-magnitude data in Figure 0) We 
compare the BSL03 models to the J-band observations, as this bandpass represents the peak of 
the near-IR SED. (Note that this is largely equivalent to using the BSL03 atmospheric models for 
their bolometric corrections.) The results are listed in Tableland agree reasonably well with those 
derived from L^oi and assuming the bolometric corrections from the known T9-T10 dwarfs. 

Measuring the differential deuterium abunda nce would help const rain the age of the system, 
analogous to the binary lithium test proposed by Liu h Leggett (j2005l ). The deuterium depletion 
timescales are long for such low-mass objects, so the exact age at which th ey become fully dep leted 
is fuzzy. Using a finely interpolated version of the Lyon/COND models (jBaraffe et al.ll2003l ) and 
constrained by the observed luminosity of each component, we find that the full range of 99% to 
0% deuterium content for the primary is 0.5-1.3 Gyr, while the the 99-0% deuterium range for 
the secondary is 2.3-5.4 Gyr. Thus, for ages of <0.5 Gyr, both components retain their deuterium 
and would have masses of <9 and <4 Mj up . Between 0.5-2 Gyr, the primary is depleted but the 
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secondary retains its deut erium. And for age s older than 2 Gyr both components are depleted and 
are >20 and >10 Mj up . (jAllers et al,ll2010l report the discovery of a young field L dwarf binary 
whose estimated luminosities and age indicate that the two components may be in the process of 
differenti al deuterium depletion .) Note that the deuterium abundance is expected to be difficult to 
measure (jChabrier et al,ll2000al ). An additional complication is the theoretical predictions for the 
D-burning limit depend on the metallicity, init ial D and He abund ances, and atmospheric boundary 
conditions and can range from sal 1-14 Mj up (jSpiegel et al.ll2010l ). 



3.5. Orbital Period 



To estimate the orbital period of the binary, we ad opt the statist i cal co nversion factor be- 
tween projected separation and true semi-major axis from lDupuv &: Liul (|201ll ). They offer several 
choices, based on the underlying eccentricity distribution and degree of completeness to finding bi- 
naries by imaging ( "discovery bias" ) . We adopt the eccentricity distribution from their compilation 
of ultracool visual binaries with high quality orbits and assume their case of modest discovery bias, 
namely that the smallest possible projected separation for detection is half the semi-major axis. 
This gives a multiplicative correction factor of 1. 07^027 (68.3% confidence limits) for converting 
the projected separation into semi-major axis (compared to 1.10^*36 f° r a uniform eccentricity 
distribution with no discovery bias). The resulting orbital period estimates range from 20-35 years 
for ages of 1.0-5.0 Gyr, albeit with significant uncertainties (Table[3j). Dynamical mass es for ultra- 



cool b inaries may be realized from orbital monitoring covering >30% of the period (e.g. iBouv et al 



2004 : iLiu et al.1 12008 : IDupuv et al.1 12009 ). Thus, CFBDSIR J1458+1013AB warrants continued 
high angular resolution imaging to determine its visual orbit, which can then be combined with the 
parallax to directly measure its total mass. 



Discussion: Y or Y not? 



Definition of the Y spectral type awaits discovery of objects whose spectra are unambiguously 
different from the T dwarfs, in the same fashion as the distinction between M and L dwarfs and 
between L and T dwarfs. The near-IR spectral progression observed for the T8-T10 dwarfs is 
subtle, as the CH4 bands are nearly fully absorbed. Such late-type objects are dist inguished from 



each other by only modest ch a nges in the shape of the near-IR flux peaks (e.g. iDelorme et al 



2008a 



Burningham et al 



more distinctive (jLeggett et al 



2008 



2010b! ). 



Lucas et al.ll2010i ). though the thermal-IR (A > 3 /im) colors are 



Recent work has highlighted near-IR spect ral features that might arise from NH3 (jSaumon et al 



2007 



Delorme et al 



2008a 



Lucas et al.ll2010l ) but not in a definitive fashion. Mod el atmospheres 



suggest the onset of NH3 a t «600 K but also indicate the effect could be indistinctive (jBurrows et al 



2003 



Leggett et al.1 120071 ). The emergence of narrow intra-band features, possibly linked to NH3 
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absorption, is difficult to disentangle from the effects of differing gravities and metallicities. Also, 
the available opacities for CH4 and H2O are known to be incomplete, complicating the identifica- 
tion of new species. In short, the current census does not seem to compel the use of the Y spectral 
type. (Of course, future identification of unambiguous examples may lead to re-classification of 
some known objects as early Y dwarfs.) 

Overall, CFBDSIR J1458+1013B is the most promising candidate to date for the hypothe- 
sized Y spectral class, given its extremely low luminosity, atypical near-IR colors, and cold inferred 
temperature. While spectroscopy is needed to place it in full context with the known T dwarfs, 
the singular nature of this object is unambiguous, thanks to our parallax measurement. CFBD- 
SIR J1458+1013B is ~1. 5-2.0 magnitudes less luminous in the near-IR than any previously known 
brown dwarf and therefore is cooler than any known brown dwarflf] For reference, two m agnitudes 



in absolute magnitude spa ns most of the T dwarf sequence, from about TO to T6 (e.g. iLiu et al 



20061 : iMarocco et al.ll2010i ). On this basis alone, this object is a significant leap. For comparison, 
the T10 dwarf UGPS J0722-05 is only re0.2 mag fainter than CFBDS J0059-01 in absolute mag- 
nitude, and among the small sample of >T8 dwarfs, the differences in absolute magnitudes are 
<0.5 mag between objects. For this reason, a spectral type estimate of >T10 seems appropriate 
for CFBDSIR J1458+1013B, as opposed to >T10. (A naive extrapolation of the linear trend in 
absolute magnitude with T subclass would make the object around T12.) While it could be that 
near-IR spectra changes so slowly with temperature that a broad luminosity range of objects all 
appear to be the same as the UGPS J0722— 05, this seems unlikely given the behavior of the hotter 
T subclasses. 

Regardless of its final spectral classification, CFBDSIR J1458+1013B provid es a gateway into 



new processes in low-temperature atmospheres. Models by iBurrows et al.l (|2003l ) predict the dis- 
appearance of the very broad sodium (0.59 fim) and potassium (0.77 fim) resonance lines at 
T eff re 450 - 500 K, co mparable to the estimated temperature of CFBDSIR J1458+1013B. (See 



also lLeggett et alj|2007l .) This change might be detected from resolved far-red optical spectroscopy 
and photometry from HST. Also, Burrows et al. predict the formation of photospheric water clouds 
at ~400-500 K (with hotter temperatures for higher surface gravities), though this phenomenon 
has little impact on their model emergent spectra. At re300-400 K, they predict the cessation 
of the blueward trend in near-IR colors with decreasing temperature. Such signatures might be 
delineated with resolved multi-band far-red optical and near-IR data of CFBDSIR J1458+1013AB. 
Our current measurements indicate that the trends of near-IR colors with spectral type continue 
without enormous changes, but that they may be non-monotonic. Indeed, the subtle changes seen 
in the near-IR SEDs of late-T dwarfs, combined with the strong alkali line depletion predicted at 



From a deep Spitzer extragalactic survey, lEisenhardt et all ^Oiol ) have identified a very faint, iow/no prop 



XT 



motion source whose exceptionaily red thermal-fR coiors and non-detection in the near-IR suggest a <450 K object at 
>70 pc. Direct comparison with CFBDSIR J1458+1013B is not possible, since there are no photometric detections 
in common between the two objects. Their object's H-band non-detection (H > 24.2 mag) is consistent with 
CFBDSIR J1458+1013B, which would be H = 24.9 mag at 70 pc. 
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<500 K, suggest that far-red optical wavelengths may become valuable for spectral classification of 
the coolest new discoveries. 

The fact that CFBDSIR J1458+1013B resides in a binary system of known distance offers 
additional future opportunties for understanding the coldest brown dwarfs. The estimated mass 
ratio (q ~ 0.5) is uncommonly small compared to most field systems, representing only <5% of 
known binaries (|Liu et al.ll2010l ). Measuring accurate Lbol an d T e R for the two components could 
test the joint accu r acy of evolutionary and atmospheric models under the constraint of coevality (e.g. 
White et al.l 11999 : iKraus k, Hillenbrand! 120091 ; iLiu et al.ll2010l ). Furthermore, given their identical 
composition and similar surface gravities, comparing the SEDs of the two components directly 
probes the influence of decreasing effective temperature, without other physical parameters to 
consider. So the binary may yield valuable insights about the transition from T dwarfs to Y dwarfs. 
Finally, the estimated orbital period of ~20-30 years means a direct dynamical mass estimate is 
feasible over the next several years; such a mass benchmark would provide a stringent test of 
substellar models. A dynamical mass would also help to demarcate the low-mass limit for star 
formation. 
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A. Spectroscopic Distance Estimates for Late-T Dwarfs 

Using the fitting results from § 13.2.31 we can test the accuracy of atmospheric models by 
co mparing the measur ed distances to model-derived spectroscopic distances following the method 
of Bowler et al. (j2009l ). The scaling factor between the model emergent spectra and the observed 
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spectra is equal to (R/d) 2 , so if a radius is known (or assumed) then the distance can be computed. 
Bowler et al. estimated distances to four late-M and L subdwarfs using this method and found 
agreement to within ~10% of published parallaxes. Several subseque nt studies have used this 



technique to estimate distances for late-T dwarfs without parallaxes (e.g. iMainzer et al.ll201ll ). but 
this method has yet to be tested for such objects. We do so now with the Ames-Cond, BSL03, and 
BT-Settl-2009 models. 

To compute model-derived spectroscopic distances, for the Ames-Cond and BT-Settl-2009 
atmospheric models we use radii from the companion grid of Lyon/Cond evolutionary models 



(jBaraffe et al.l 120031 ) interpolated across ages, effective temperatures, and surface gravities. For 
the BSL03 models , we u se a power law fit for the radius as a function of T c g and log (g) from 
the Burrows et al. (jl997l ) evolutionary models (Equation 1 of BSL03). We compute errors in the 



spectroscopic distances from two components. (1) Measurement uncertainties from the finite S/N 
of the spectra and the absolute flux calibration from the J-band photometry are accounted for 
in a Monte Carlo fashion. We adopt the median and RMS of the resulting distributions as one 
version of the results. (2) In addition, we account for systematic uncertainties in the models by 
assuming normally distributed errors of 50 K and 0.25 dex in the fitted T e g and log(g). The 
standard deviation of the resulting spectroscopic distance distribution is reported along with the 
uncertainty due to random errors. 

Table [5] summarizes the results, with both sets of uncertainty calculations. The agreement 
between the spectroscopic distances and the parallactic distance range from «10% to about a factor 
of two, with no apparent trend in accuracy with distance, spectral type, or model grid. The large 
scatter likely arises from the sensitivity of the distance estimate to the best-fitting model T e g. The 
typically larger spectroscopic distances compared to parallactic distances imply an over-estimate of 
T e fj in the best-fitting model, because the scaling factor is highly sensitive to temperature (absolute 
flux level) and very weakly dependent on gravity. 
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Fig. 1 — Keck LGS AO imaging of CFBDSIR J1458+1013AB. North is up and east is left. Each 
image is 0.75" on a side. The greyscale images are shown with a power-law stretch with an exponent 
of 1/4. The contours are shown after smoothing the images by a 2-pixel gaussian. Contours are 
drawn from 90%, 45%, 22.5%, 11.2%, and 5.6% of the peak value in each bandpass. 
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Fig. 2.— The change in the position of CFBDSIR J1458+1013B relative to 

CFBDSIR J1458+1013A. The red point shows the measured change from our two Keck 
LGS epochs, which is consistent with no motion. The blue symbol shows the expected motion (and 
its uncertainties) if component B is a background source, given the measured parallax and proper 
motion of CFBDSIR J1458+1013. (Over the 1.5 month time baseline between Keck epochs, the 
proper motion dominates.) The two points are inconsistent at 9.1a, indicating that the system is 
a physical binary. 
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Fig. 3. — Left: Relative astrometry from CFHT/WIRCam J-band imaging of 
CFBDSIR J1458+1013AB that shows both proper motion and parallax. The origin corre- 
sponds to the first epoch of data. Middle and Right: Top panels show relative astrometry in a and 
5, respectively, as a function of Julian year after subtracting the best-fit proper motion. (This is for 
display purposes only; in our analysis we fit for both proper motion and parallax simultaneously.) 
Bottom panels show the residuals after subtracting both parallax and proper motion. 



-23- 



1.0- 



0.5- 



"D 

CD 

J 0.0 
E 




1.0 



— i — i — i — i — i — 

UGPS0722 (T10) 
CFBDS1458AB (T9.5) 

ULAS0034 (T9)- 
ULAS1335 (T9) 
CFBDS0059 (T9) 

CIA H, 



CH 4 



H ? 



1.5 



2.0 



2.5 





1.0 1.1 1.2 1.3 



1.5 



1.6 



1.7 



X (jim) 



Fig. 4 — The integrated-light near-IR spectrum of CFBDSIR J1458+1013AB (black line) com- 
pared with other late-T dwarfs. The spectra have been smooth ed for clarity. The upper panel shows 
a comparison with the T10 dwarf UGPS J0722— 05 (red line; iLucas et al.ll2010l ). The lower pan- 
els show the Y J (left) and iJ-band (right) comparison with UGPS J0722— 05 and three T9 dwarfs 



(Warren et al 



2007 



Burningham et al. 



2008 



Delorme et al.ll2008al ) . The black horizontal bars near 



the bottom of the plots show the wavelength ranges used to normalize the spectra (1.25-1.28 /jm 
for the top and bottom left plots, and 1.57-1.59 fim for the bottom right plot). Note that for the 
upper panel, the x-axis uses a log scale to help make the bluer bandpasses more visible. 
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Fig. 5. — The Wj and J/K spectral indices for very late-T dwarfs compared to the solar-metallicity 
BT-Settl-2009 model atmospheres. (See Table 0] for measurements and references.) The model grid 
has b een scaled to agree with the T8.5 benchmark dwarf Wolf 940B (T e g = 575 K, log(g) = 
4.75; iBurningham et al.l |2009i ). The plotted errors represent the la uncertainties in the index 
measurements. 
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Ames-Cond: 800.0/4.5 
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Fig. 6. — Best-fitting BT-Settl-2009 (blue), Ames-Cond (red), and BSL03 (green) model atmospheres to 
late-T dwarfs with parallaxes. (See also Table [5]) The gray horizontal bars show the wavelengths used in 
the fit, chosen to avoid known systematic errors in the models (e.g., the Y"-band and the red side of the H- 
band) and regions of low S/N. The observed spectra are normalized to unity flux between 1.25-1.28 /j,m, and 
the normalization of the models is an outcome of the fitting process. The observed spectra for Ross 458C, 
Wolf 940B, CFBDS J0059-01, ULAS J0034-00, and ULAS J1335+11 arc smoothed for clarity here, but 
the models were fit to the original (unsmoothed) data. Note that the x-axis uses a log scale, to help make 
the bluer bandpasses more visible. 
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Fig. 7. — Color-magnitude diagrams comparing CFBDSIR J1458+1013AB (red dots) to field L and T dwarfs with parallaxes 



(jDahn et alJl2002l ; iTinney et al 



2003; 



Vrba et al. 



2004: 



Marocco et al.ll2010l ). based on the compilation by lLeggett et al.1 (|2010bl ) 



and plotted a s filled diamonds. Known binaries are exclude d. T he TIP dwarf UGPS J0722— 05 is labeled. We also show the 
Lyon/COND (jBaraffe et al.ll2003h . ISaumon & Marlevl fcoosh . and burrows et all (|2003h models. For the SM08 models, we plot 
the log(^) = 5 models that have a moderate cloud layer (f sec i = 4) and those that have no clouds ("nc"). The open circles 
on the model tracks demarcate a common set of T e g 's; for models plotted with a solid line, the corresponding T c g values are 
listed on the right hand axis for the Saumon & Marley models (which go down to 500 K) and then the BSL03 models for cooler 
temperatures. This combination is chosen for the labels because it agrees better with the data than the Lyon/COND models, 
though in general no single model agrees well with the loci of the latest T dwarfs and CFBDSIR J1458+1013B. 
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Fig. 8. — Comparison of the colors and magnitudes of CFBDSIR J1458+1013 A and B to those 
of field objects. The latter come from the compil ation of iLeggett et al.l (l2010bh. augmented with 



resolved photo metry of tight ultracool binaries (]McCaughrean et al.l 12004 ; iLiu Sz Leggett 
Liu et alJbooi ). 
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Fig. 9. — Compa rison of model-derived spectroscopic distances (d spec ) using the method of 
Bowler et alJ (I2009I ) with the parallactic distances (d n ) for the sample of late-T dwarfs in Table EJ 
See Appendix for details. 



Table 1. Keck LGS AO Observations 



Date 


Filter 3 - 


Airmass 


FWHM 


Separation 


Position angle 


Amag 


(UT) 






(mas) 


(mas) 


(deg) 




2010-May-22 


CHiS 


1.14 


62 ±6 


108 ±6 


308 ±2 


2.17±0.12 


2010-Jul-08 


J 


1.33 


94 ± 10 


110 ±11 


308 ±4 


1.8 ±0.4 




H 


1.04 


61 ±4 


111 ±3 


309.3 ± 1.4 


2.33 ±0.14 




K 


1.11 


60 ±3 


122 ±6 


314 ±2 


2.2 ±0.2 



Note. — All photometry on the MKO system. The tabulated uncertainties on the imaging 
performance (FWHM and Strehl ratio) are the RMS of the measurements from the individual 
images. 



Table 2. CFHT/WIRCAM J-band Astrometry of CFBDSIR J1458+1013 



N cp N lei Epoch 1 


At (yr) 


AAM 


N {l (N cal ) 


\x (mas/yr) 


P. A. (°) 7T re i (mas) 


Air (mas) 7r a b s (mas) 


X 2 (dof) 


7 72 2009.42 


1.07 


0.03 


68 (8) 


432 ±6 


154.2 ±0.7 42.4 ±4.5 


0.9 ±0.1 43.3 ±4.5 


11.5 (9) 



Note. — N ep = Number of distinct observing epochs (i.e., nights); 7V rc f = Number of reference stars used; Epoch 1 = Date of our first 
WIRCam observation; At = Time baseline for the observations; AAM = Maximum difference in airmass over all epochs; Nf T = Total number 
of frames obtained (10-11 per epoch); A^ ca i = Subset of reference stars used in the absolute calibration; Air = Correction applied to the 
best- fit parallax to determine the absolute parallax (7r a b s = 7r re i + An) . 
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Table 3. Properties of CFBDSIR J1458+1013AB 



Property 



Component A 



Component B 



Integrated-light spectral type 
d (pc) 

Projected separation (AU) 
Vtan (km/s) 



T9.5 
23.1 ± 2.4 
2.6 ± 0.3 
47 ± 5 



J (mag) 


19.86 ± 0.07 


21.66 ± 0.34 


CH4S (mag) 


19.73 ± 0.10 


21.90 ± 0.15 


H (mag) 


20.18 ± 0.10 


22.51 ± 0.16 


K (mag) 


20.63 ± 0.24 


22.83 ± 0.30 


J — H (mag) 


-0.32 ± 0.12 


-0.85 ± 0.37 


J — K (mag) 


-0.78 ± 0.25 


-1.18 ± 0.45 


CH 4 s - H (mag) 


-0.45 ± 0.14 


-0.61 ± 0.22 


H - K (mag) 


-0.45 ± 0.26 


-0.32 ± 0.34 


M(J) (mag) 


18.04 ± 0.23 


19.84 ± 0.40 


M(CH 4 s) (mag) 


17.91 ± 0.24 


20.08 ± 0.27 


M(H) (mag) 


18.36 ± 0.24 


20.69 ± 0.27 


M{K) (mag) 


18.82 ± 0.33 


21.02 ± 0.37 


Near-IR spectral type 


T9.5 


>T10 


log(L bol /L Q ) (dex) 


-6.02 ±0.14 


-6.74 ±0.19 



Evolutionary model results for 1.0 and 5.0 Gyr: Lyon/COND models and Lt 



Mass {M Jup ) 
Total mass (Mj up ) 
q (ee M B /M A ) 

l°g(g) (cgs) 
Orbital period (yr) 


12.1 ± 1.9, 31 ± 4 5.8 ± 1.3, 14 ± 3 
18 ± 2, 45 ± 5 
0.48 ± 0.14, 0.45 ± 0.12 
556 ± 48, 605 ± 55 360 ± 40, 380 ± 50 
4.45 ± 0.07, 5.00 ± 0.08 4.10 ± 0.10, 4.58 ± 0.11 

OC- + 28 nri-^-lS 

OJ -10> ^ -6 


Evolutionary model results for 1.0 and 5.0 Gyr: Burrows et al. (1997) models and Lbol 


Mass (M Jup ) 
Total mass (Mj up ) 
q (ee Mb/Ma) 

Teff(K) 

log(9) (cgs) 
Orbital period (yr) 


13 ± 2, 36 ± 4 6.8 ± 1.5, 17 ± 4 
20 ± 3, 53 ± 6 
0.53 ± 0.15, 0.47 ± 0.13 
550 ± 50, 600 ± 60 350 ± 40, 380 ± 50 
4.47 ± 0.07, 5.06 ± 0.07 4.14 ± 0.10, 4.65 ± 0.12 
33± 9 7 , 20lg 7 


Evolutionary model results for 1.0 and 5.0 Gyr: Burrows et al. (2003) models and M(J) 


Mass {M Jup ) 
Total mass (Mj up ) 
q (ee Mb 1 Ma ) 

log(ff) (cgs) 
Orbital period (yr) 


11.1 ± 0.7, >25 7.6 ± 0.6, 18.8 ± 1.3 
18.7 ± 0.9, >44 
0.68 ± 0.07, <0.75 
479 ± 20, >483 386 ± 15, 407 ± 15 
4.37 ± 0.03, >4.85 4.19 ± 0.04, 4.69 ± 0.03 
34^0, <22 
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Note, — A l l infr ared photometry on the MKO photometric system. For the 5 Gyr 
iBurrows et al. I <2003h evolutionary model fits, component A is brighter in M(J) than the 
model grid, so its derived physical parameters are listed as limits. 



Table 4. Spectral indices for very late-T dwarfs. 



Object 


Sp. Type 


H 2 0-J 


Wj 


CH 4 -J 


H 2 O-H 


CH 4 -H 


NH3-H 


CH4-K 


K/J 


Wolf 940B 


T8.5 


0.030 ± 0.005 


0.274 ± 0.002 


0.1521 ± 0.0009 


0.142 ± 0.002 


0.090 ± 0.002 


0.537 ± 0.002 


0.072 ± 0.013 


0.1113 ± 0.0014 


ULAS J0034-00 


T8.5 


0.012 ± 0.006 


0.262 ± 0.004 


0.145 ± 0.009 


0.133 ± 0.010 


0.097 ± 0.006 


0.510 ± 0.011 


0.092 ± 0.015 


0.128 ± 0.003 


CFBDS J0059-01 


T8.5 


0.029 ± 0.005 


0.257 ± 0.004 


0.165 ± 0.005 


0.119 ± 0.008 


0.084 ± 0.002 


0.526 ± 0.005 


0.128 ± 0.037 


0.101 ± 0.002 


CFBDSIR J1458+1013AB 


T9.5 


-0.002 ± 0.008 


0.237 ± 0.006 


0.103 ± 0.009 


0.078 ± 0.016 


0.044 ± 0.018 


0.54 ± 0.02 


0.06 ± 0.19 


0.058 ± 0.012 


UGPS J0722-05 


T10 


0.0339 ± 0.0017 


0.2074 ± 0.0012 


0.1358 ± 0.0019 


0.1218 ± 0.0017 


0.0643 ± 0.0013 


0.492 ± 0.002 


0.096 ± 0.002 


0.1262 ± 0.0003 


Note. — Near-IR. spectral 


indices and 


provisional spectral types for some of the 


latest-type brown 


dwarfs SWarren et al. 


2007: Burninsham et al. 2009: Delorme ct al. 2008s 


hlLucas et alj|201di. 


The index definitions are from iBureasser 


et al.l l2006bH.IWarren 


et al.l J2007I1. andl 


Delorme et al. (2008ah. Note that the ni 


jgativc value of H2 


O-J for CFBDSIR J1458+1013 is 


due to the fact that 


there is no detected flux in the 1.14—1.165 fim water absorption 


band probed by this index. 













I 

CO 
CO 

I 



Table 5. Atmospheric Model Fits to Late-T Dwarfs with Parallaxes 



Amcs-Cond BSL03 BT-Settl-2009 









log 9 


iT 


<^spcc (fincj "^tot ) 


^eff 


log 9 


BP 


<^spcc (°~mc 


\D 

O"tot ) 




log 9 


R c 


^spec 


7 \D 

(cr mc , (T t ot ) 






Name 


Ref 


(K) 


(cgs) 


(Rjup) 


(PC) 


(K) 


(cgs) 


(Rjup) 


(pc) 




(K) 


(cgs) 


(Rjup) 




(PC) 


(PC) 


7T Rcf 


Gl 570D (T7.5) 


1 


800 


4.5 


1.06 


6.8 (0.1, 1.4) 


797 d 


4.80 d 


0.996 


14.8 (0.2, 


2.6) 


900 


5.0 


0.903 


7.9 


(0.2, 1.4) 


5.84 (0.03) e 


13 


HD 3651B (T7.5) 


2, 3 


800 


1.5 


1.06 


13.1 (0.2, 2.5) 


797 d 


4.80 d 


0.996 


14.1 (0.2, 


2.1) 


850 


1.5 


1.07 


15.8 


(0.2, 2.3) 


11.06 (0.04) 


13 


2MASS J0415-09 (T8) 


4 


800 


5.0 


0.89 


6.29 (0.09, 2.4) 


703 


4.07 


1.33 


6.50 (0.1, 


2.9) 


600 


5.5 


f 




f 


5.65 (0.09) 


11 


Ross 458C (T8.5) 


5, 6 


900 


4.0 


1.24 


23.1 (0.09, 3.2) 


859 d 


4.24 d 


1.27 


24.2 (0.1, 


2.5) 


850 


4.5 


1.07 


17.8 


(0.1, 2.6) 


11.7 (0.2) 


13 


Wolf 940B (T8.5) 


7 


800 


1.5 


1.06 


30.4 (0.3, 5.5) 


483 d 


4.85 d 


0.924 


8.37 (0.07 


2.9) 


650 


5.0 


0.874 


15.5 


(0.1, 3.7) 


12.5 (0.7) 


15 


CFBDS J0059-01 (T9) 


8 


700 


4.5 


1.04 


20.2 (0.3, 4.1) 


483 d 


4.85 d 


0.924 


8.01 (0.1, 


2.7) 


600 


5.0 


0.879 


11.5 


(0.2, 3.3) 


9.2 (0.4) 


16 


ULAS J0034-00 (T9) 


9 


800 


4.5 


1.06 


30.5 (0.5, 5.8) 


686 


4.66 


1.04 


24.7 (0.4, 


5.8) 


750 d 


4.0 d 


1.17 


30.7 


(0.5, 4.5) 


12.8 (0.6) 


16, 17 


ULAS J1335+H (T9) 


10 


800 


1.5 


1.06 


26.4 (0.2, 5.2) 


686 


4.66 


1.04 


20.6 (0.1, 


5.0) 


800 d 


4.0 d 


1.21 


29.9 


(0.1, 4.6) 


10.3 (0.3) 


16 


CFBDSIR J1458+1013AB (T9.5)« 


11 


700 


4.5 


1.04 


46.8 (1.0, 10) 


555 d 


4.83 d 


0.946 


27.2 (0.9, 


7.6) 


700 


5.0 


0.880 


40.0 


(0.8, 9.3) 


23.1 (2.4) 


18 


UGPS J0722-05 (T10) 


12 


700 


4.5 


1.04 


9.9 (0.1, 2.0) 


421 d 


4.72 d 


0.961 


2.62 (0.02 


1.4) 


600 


5.0 


0.879 


5.6 


(0.1, 1.6) 


4.1 (0.5) 


12 



Note. — Best-fitting Amcs-Cond iAllard et all200ll) . BSL03 jBurrows et all2003h . and BT-Settl-2009 iAllard fc Frevtaell201(f) model spectra to the combined 1.15-1.35 fim., 1.45-1.6 fJLta, 
and 1.95-2.25 fim regions of the observed spectra. 

a Radius derived from an interpolated grid of Lyon/Cond evolutionary models. This value corresponds to the best-fitting T g and log g Ames-Cond atmospheric model. 

b Model-derived spectroscopic distance estimate based on the method of iBowler et al. | l l2009l) . &mc denotes the Monte Carlo-derived distance error associated with the best-fitting 
atmospheric model and incorporates measurement and flux-calibration errors. <x tot is the Monte Carlo-derived total measurement error in the distance, based on the measurement errors 
and assuming lcr uncertainties of 50 K and 0.25 dex in T e ff and log g in the model atmosphere fits. 

c Radius computed using the best-fitting T e g and log g and the power law fit to R(T c g, log g) from Equation 1 of lBurrows et all £2003). 

d Bcst-fit model is at the edge of the model grid. 

°Distance to the primary Gl 570 A. 

f The best-fitting model atmosphere parameters are outside the grid of evolutionary models, so there is no corresponding radius or distance estimate, 
integrated light spectrum. 



Ref eren ces. — l-|Burgasser et all ibOOOT); 2-|Mugrauer et TaLl 1120061) ; 3-|Luhman et al. J2007h; 4-|Burgasser et all j2002|); 5-lGoldman etakl bold); 6-|Scholz 1 20101: 
<2009l); 8-|Delorme etal. j20q8ah : 9-|Warren et al. (120071) ; 10— Burningham et al. I <200a) ; 11- lDelorme et alj J20ld ); 12- |Lucas et al, I j20ld) ; 13- lvan Leeuwenl J2007ll 
15- |Harrington fc Dahnl il980h 16- |Marocco et alj J20ldfh 17- ISmart et al.l ||2010|~ 18-this work. 



-|Bjmnn^lianie^alJ 
4Vrba et al. I j2004h 



-35- 



Table 6. Physical Properties Derived from Model Atmospheres and Parallactic Distances 
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log 9 
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lot 


! 9 
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Name 


(K) 


(cgs) 


(M Jup ) 


(K) 


(cgs) 


(M Jup ) 


(K) 


(cgs) 


(M. 


lup) 


Gl 570D (T7.5) 


780- 


-880 


4.7- 


-5.4 


21 


51 


>520 


>4.4 


>13 


760- 


-880 


4.7- 


-5.4 


21 


54 


HD 3651B (T7.5) 


780- 


-880 


4.7- 


-5.4 


21 


54 


>700 


>4.6 


>22 


740- 


-860 


4.7- 


5,1 


20 


54 


2MASS J0415-09 (T8) 


790- 


-800 


4.7- 


-5.3 


23 


53 


>660 


>4.6 


>19 


600- 


-720 


4.5- 


5.2 


9 


45 


Ross 458C (T8.5) 


690- 


-760 


4.0- 


-4.7 


6 


20 


>590 


>3.9 


>5 


660- 


-740 


4.0- 


-4.7 


6 


20 


Wolf 940B (T8.5) 


570- 


-680 


4.5- 


-5.2 


12 


43 


>500 


>4.4 


>12 


560- 


-640 


4.4- 


5.2 


12 


12 


CFBDS J0059-01 (T9) 


520- 


-600 


4.4- 


-5.2 


10 


35 


>440 


>4.3 


>10 


500- 


-600 


4.4- 


5.1 


11 


34 


ULAS J0034-00 (T9) 


560- 


-680 


4.5- 


-5.2 


12 


13 


>500 


>4.4 


>12 


560- 


-660 


4.4- 


5.2 


12 


42 


ULAS J1335+11 (T9) 


560- 


-660 


4.4- 


-5.2 


12 


42 


>500 


>4.4 


>12 


560- 


-640 


4.4- 


5.2 


12 


12 


CFBDSIR J1458+1013AB (T9.5) a 


520- 


-640 


4.4- 


-5.2 


11 


42 


>460 


>4.3 


>10 


520- 


-620 


4.4- 


5.2 


11 


36 


UGPS J0722-05 (T10) 


480- 


-600 


4.4- 


5.1 


10 


35 


>420 


>4.3 


>9 


480- 


-580 


4.4- 


-5.1 


10 


33 



integrated light spectrum. 



